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Abstract 

The charge transfer from an adatom to a semiconductor substrate of one-dimensional quantum 
dot array is evaluated theoretically. Due to the Van Hove singularity in the density of electron 
states at the band edges, the charge transfer decay rate is enhanced nonanalytically in terms of 
the coupling constant g as g 4//3 . The optical absorption spectrum for the ionization of a core level 
electron of the adatom to the conduction band is also calculated. The reversible non-Markovian 
process and irreversible Markovian process in the time evolution of the adatom localized state 
manifest themselves in the absorption spectrum through the branch point and pole contributions, 
respectively. 
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I. INTRODUCTION 



Thanks to recent advances in nanotechnology,various types of artificial low-dimensional 
semiconductor structures have been fabricated The quantum confinement of electron 
in these structures greatly modifies the density of states of carriers resulting in complete 
different electronic and optical properties from the bulk system . Recently one dimen- 
sional quantum wire and quantum-dot array have been manufactured in various ways 
and the formation of one dimensional miniband has been theoretically and experimentally 
investigated [a El Q • It has been found that the Van Hove singularity of the density of state 
inherent with the one-dimensionality causes a characteristic electronic transport 3 0] . 

In this paper, we consider the charge transfer between an adatom localized state and 
a one dimensional miniband associated with a quantum wire or quantum-dot array. The 
charge transfer between an adatom and a substrate semiconductor has been extensively 
studied [10] . We will show that the low dimensionality greatly modifies the charge transfer 
process from the adatom to the semiconductor quantum dot array due to the singularity of 
the density of states. The physical situation we consider in this paper is shown in Fig. da) 
where an adatom is attached to a semiconductor quantum-dot array surface (hereafter we 
simply call it a superlattice). 

The charge transfer of an electron from adatom to the miniband of the superlattice is 
caused by the hybridization of the adatom wave function with the miniband which can 
be described by the bilinear coupling between the adatom localized state and the bound 
state of a single quantum dot in which the adatom is situated. The situation may be 
described by one-dimensional version of Newns- Anderson model which has been extensively 
used to investigate the charge transfer process between the adatom and the semiconductor 



substrate 



[lO, 



111 ] . As shown below, the one-dimensional Newns Anderson model we consider 



here is equivalent to the Friedrichs Hamiltonian which we have investigated in our previous 
letter [l2|. 

In this recent letter we reported a vast increase in the decay rate of an excited dipole 
molecule traveling in an one-dimensionally confined electromagnetic waveguide when the 
cutoff frequency of the waveguide is near the characteristic frequency of the dipole jl^ . 
This vast increase is a direct consequence of a singularity in the density of photon states at 
the cutoff frequency. Due to this singularity, standard perturbation analysis breaks down 
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and hence one cannot apply Fermi's golden rule to evaluate the decay rate in the vicinity of 
the singularity We have shown that in this case the decay rate of the excited dipole depends 
non-analytically as g 4//3 on the coupling constant g. In the present article, we report that 
the same nonanalytic enhancement of the decay rate (including the g 4 / 3 law) can be found 
in our system, despite the fact that the density of electron states of the miniband of the 
superlattice is different from the density of states of photon states in a waveguide. 

Although the exponential decay law for the unstable state has been observed ubiquitously 
in nature, quantum mechanics predicts that there should be a deviation from the exponential 
decay law for the unstable state 3|. It has been shown that, irrespective of any specific form 
for the interaction, the time evolution of the surviving amplitude of an unstable state deviates 
from the exponential law on short and long time scales due to the existence of the lower 
bound on the energy, i.e., the branch point effect. Wilkinson et al. has recently succeeded 
in measuring the branch point effect in super-cooled sodium atoms even though the 
timescale in which the deviation from the exponential decay law appears is very short, Thus 
separation of the pole and branch point effects in the time evolution of the unstable state is 
essential to our understanding of the decay process, as has been done for a system composed 



of an excited atom coupled to a radiation field [la, [16( . The separation of the two effects is 
also useful because the reversible non-Markovian process in the decay process (due to the 
braM h point effeet, is dffeeti, reiated t o the q _ Z» effect ffilQQ. 

Our goal is to present an actual experimental method that enables us to separately 
evaluate the pole and branch point effects for the decay process of the unstable state in this 
system. In this article, we consider the optical absorption process in which a core electron 
of the adatom with discrete energy is transferred to the continuous conduction band. It will 
be shown here that the spectral shape is influenced by the Markovian process due to the 
pole effect and also by the non-Markovian process due to the branch point effect. 

The Friedrichs model presented here is equivalent to that known as the Fano model which 
was originally developed to explain the absorption spectrum for the autoionization process 
in the He atom j^]. It is well known that this absorption spectrum has an asymmetric 
spectral profiles due to quantum interference between different optical transition paths. 
The appearance of the quantum interference indicates that the quantum coherence, which is 
a source of the memory effect, plays a key role in the decaying process of the excited state. 
As will be shown in the present paper, since the branch point effect which accounts for the 
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non-Markovian decaying process with a memory effect is intensified by the singularity in 
the density of states, the absorption spectral profile provides us with the information on the 
extent of the persistence of the quantum coherence in the decaying process. 

In §2 we present the model and show the nonanalytic enhancement of the decay rate for 
the unstable state. We decompose the absorption spectrum into contributions from the pole 
and branch point effects in §3. In §4 we summarize our results and provide some discussion. 



II. MODEL AND NONANALYTIC ENHANCEMENT OF DECAY RATE 

We consider a ID semiconductor superlattice with an adatom on the surface, as shown 
in FiefTTa). The width of each potential well is a few nm to 100 nm. One can make this 

n n 

device, for example, with GaAs/GaAlAs heterostructures [1|,|5|. The superlattice consists of 
N ^> 1 identical quantum wells, and each well is assumed to have a single bound state |n) 
(n = — N to N) of equal energy, where the tilde is used to distinguish the site representation 
from the wavenumber representation below. A miniband is formed in the superlattice due 
to the electron tunneling through the potential barrier j^, 0, Q, 13 • We assume only nearest 
neighbor tunneling occurs with a transition probability of —B/2. The ID superlattice is 
then represented by the one- dimensional tight binding model, and we have a continuous 
miniband of width 2B in the limit of iV — > oo. In addition to the miniband we consider 
an adatom localized state \d) with energy Eq, also we consider an inner core level |c) with 
energy E c . Both of these states are associated with the adatom impurity located at the 
n = site. The adatom localized state \d) is hybridized with the |0) state with coupling 
strength gB. 

Taking H — 1 hereafter, the electronic Hamiltonian He is then written as 
H E = E c \c)(c\ + E \d)(d\- - ^ \m)(n\ 

<m,n> 

+gB{\d)(b\ + \0)(d\) , (1) 

where < • ■ ■ > means taking nearest neighbor sum in Eq.l JIj) . By introducing the wavenumber 
representation 

1 N 

\ k ) = -/^^^[iknWn) , (2) 



n=l 
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we can rewrite He in the form of the Friedrichs model 

H E = E c \c)(c\ + E d \d)(d\ + J2 E k\ k )( k \ 



)=52gB(\d)(k\ + \k)(d\), (3) 



N k 

where E k = —Bcosk. We impose periodic boundary conditions, leading to k — 2nj/N, 
where j is an integer running from —N/2 to N/2. The energy dispersion relation gives a 
divergence in the density of states at either band edge: 

= \~w=Wi ' (4) 

For the ID Friedrichs model, the electronic Hamiltonian He is diagonalized in terms of 
the so-called Friedrichs solution which results in the spectral representation |2l[ 

2 

H E = E c \c)(c\ +J2Ei\<f>i)(<f>i\ +J2 E k\<f>k)(<l>k\ , (5) 
i=i k 

where (i = 1,2) are two stable eigenstates with energies E iy and \(p k ) with energy E k 
(where \E k \ < B) is given by 

I0D = l*> + 1 <)B 



Nv + (E k 



\\d) + -^Y — -\k')\ , (6) 



with 



= (d\ ,_ TT ,J d) 



r]+(z) ddX ' v \z-H 



p lv- g 2 B 2 



(?) 



and a positive infinitesimal e. The Green function G dd (z) is analytically continued from 
the upper half complex Ek plane to the lower half plane. Here, we do not write explicit 



forms of \4>i), since we will not use them in this paper 21]. In the limit of N — > oo, the 
summation over the wavenumber turns to the continuous integral, and then the self-energy 
term in Eq.((7J) reads 

g 2 B 2 1 r „ g 2 B 2 



1 E 



dk- 



N z-E k 2tt J_ n z-E k 

d E k p{E k )-^- = 9 . (8) 

b z-E k Vz 2 - B 2 1 ; 



By substituting the right-hand side of Eq.flHJ) into r] + (z) and taking the square of the 
dispersion relation we obtain a quartic equation in z. From the form of the quartic equation, 
it is readily seen that the solutions are symmetric about the origin E = 0. By applying 
the standard method for solving a quartic equation, one can find the explicit solutions of 
the dispersion equation r] + (z) = 0. The solutions consist of the pole Zd = E — 27 of G^ d (z) 
in the lower half complex plane in the second Riemann sheet corresponding to the unstable 
decaying state, and the poles on the real axis corresponding to the stable states, Ei{i = 1,2). 
For illustration, we plot in Fig|2fa) E and Ei for g = 0.5 as a function of Eq/B with thick 
dashed and solid lines, respectively. For arbitrary E , there always exist two stable solutions 
outside the electronic miniband. Note that there is a critical value E 1 at which the unstable 
state with imaginary part 7 7^ appears. The unstable state exists for all \Eq\ < E^, as 
indicated by the arrows in Fig|2fa). In Fig. |2fb) we plot the decay rate j/B of the unstable 
solution for 3 = 0.1 as a function of Eq/B. The maximum value of 7 = 7 ma;c occurs at 
E = ±B. 

Since the explicit form of the full solution is complicated, we present an approximate 
calculation of the maximum value 7 max of 7 at Eq/B = ±1 and the critical value E 1 of Eq 
where the unstable solution appears. To estimate •y ma x we put E = B in the dispersion 
equation. After a simple manipulation we obtain 

4/3 

c = 1 + ( - 1)2/3 TcTW' (9) 

where ( = z/B. The zeroth order solution (for g — 0) is £ = 1. We use this as our starting 
point and solve iteratively to find 

v% 4/3 g s 3 Q ( 16/3N / 10 ) 

B 2 4 /3 28/^ h 1 ' 

where the third order contribution (g 4 ^ 3, with j = 3) vanishes. 

To estimate E 1 near Eq = B we put z = Eq + g a z\ where a > and z\ are unknown 
variables which are independent of g. Then, by keeping only the predominant contribution 
to the dispersion equation, we obtain 

C - E « ~y=^= t- - , (11) 



where E = Eq/B. Squaring this equation yields 

f 

E + l 



/(C) = (C-^o) 2 (C-i)--X^o. (12) 



The function f(() represents a cubic curve. By taking the derivative of f(() and setting 
f'{() = 0, we find the threshold value of Eq = E^/B > 1 at which the complex solutions of 
/(C) = appear. This leads to the first two terms in the following expression for E y : 

-L = l + JL-- y — + 0{g 12 / 3 ). (13) 

For precision, we have presented the second order correction that is obtained from the exact 
solution of the original quartic equation. The values of 7 and E 1 are non-analytic in g for 
g — 0, and hence one cannot obtain these results from a series expansion in g using ordinary 
perturbation analysis. 

It should be noted that the g 4 / 3 law for the decay rate is rather universal around the 
edge of the continuous spectrum. Indeed, if the integration over the wavenumber in the self- 
energy is a Cauchy integral that is a double- valued function on z, as is the case in Eq.(jZJ), 
then the edge of the miniband will give a square root type of essential singularity in the 
self-energy. In addition, if this singularity leads to a divergence at the edge, one can show 
that the g 4//3 effect appears in the vicinity of the singularity. Indeed, this is also the case for 
the dipole molecule traveling inside the waveguide mentioned above [l2|. 

In terms of the solutions of the dispersion equation obtained here, we may evaluate the 
time evolution of the surviving amplitude of the adatom localized state defined by P(t) = 
\(d\ exp[— iH E t]\d)\ 2 . Indeed, aside from the persistent oscillation attributed to the two 
stable states, there are two contributions to the time evolution of the surviving amplitude; 
one comes from the pole contribution over the wavenumber integral, where the location 
of the pole is given by the solution of the dispersion relation r) + (Ek) = in the complex 
Ek plane discussed above, and the other comes from the two branch point contributions 
that are located at the edges of the miniband at ±B. As usual, the pole contribution 
leads to an exponential decay exp[— i(E — fy)t], while the branch point leads to a power 
law decay. If the condition 1 — \E \/ B ^> g is satisfied, we have a well separated time 
scale between the exponential and power law decay in the decay process of the localized 
state. The exponential law for the decay process is thus a good approximation for a time 
scale on the order of t ~ I/7. However, once this condition is no longer fulfilled, the time 
separation becomes obscure and we must take into account the branch point effects in the 
time evolution. It is, however, a cumbersome task to evaluate the integral associated with 
the branch point in the time evolution. In the next section, we will evaluate the optical 
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absorption spectrum instead, to separate the contribution of the Markovian effect from the 
non-Markovian effect. It should be generally easier in experiment to observe the absorption 
spectrum than a time evolution of the state \d>. 

III. ABSORPTION SPECTRUM 

A. Spectral representation and the separation of the pole and branch point con- 
tributions 

Here we consider the optical absorption by a transition from an inner core level of the 
adatom to the conduction electronic states of the ID semiconductor superlattice, as shown in 
Fig^b). There are two optical transition paths from the core level |c) with large transition 
amplitudes: One is that from the core level |c) to the localized state \d) and the other is 
from the core level |c) to the h = site |0). The transition probabilities are denoted by 
Tdc and Tq c , respectively. The contributions of the transitions to the other bound states |n) 
are small, so we neglect the contributions in this paper, though it is not difficult to include 
those effects. 

We consider here a single mode of the optical light field with frequency Q. The total 
Hamiltonian is then written as 

Htotai = H R + H E + H RE , (14) 

where the electronic Hamiltonian He is given in Eq.fjHJ) and 

Hr = Qa^a , (15a) 
Hre = (T dc \d)(c\+T 0c \b)(c\)a + R.c. (15b) 
= fa + a)f^ . (15c) 

The Hr represents a monochromatic light field where a(a') is the annihilation (creation) 
operator of the field, and Hre describes the interaction between the light field and the 
electronic system under the rotating wave approximation. 

In the initial state of the absorption process, we have a core electron in the |c) state and 
a single photon with energy Q; we write this initial state as |c; 1) with energy E c + Q, where 
1 denotes the photon number for the mode Q. After the absorption process the energy of 



S 



the photon is transferred to the electronic system. The final states then take the form 0) 
(i = 1, 2) or \<f>%] 0), with the energies E{ or E k respectively. 

When the value of E c + Q falls in between —B and B, the state \c; 1) is resonant with 
the electronic continuum As a result, the energy transfer from the light field to 

the electronic system occurs. Hereafter we take — E c as the origin of the light energy: 

n + e c -> a 

The decay rate of |c; 1) is then determined by the pole location of the Green's function 
for \c; 1) as 



Gt(z) 



(c;i| 



l 



|c;l). 



(16) 



— H tota i] 

The decay rate 7 a fe s (fi) is a function of f2, and thus we identify the absorption spectrum as 
F(fl) = 7 a 5 S (fi). In the weak coupling limit of Ta c and T 0c , the absorption spectrum F(Q) 
reduces to 

\(S?Af\c) 2 

F(n) = - lim Im 



— lim Im > 



O — Ek + it 
g (T dc |cj)+T 0c |0))| 2 

The matrix element in Eqs. f|17|) may be obtained by using Eq.(JS]): 

,n,F\ 9 B 1 
W fc ) = 7^ZT73 

(o|0D = i + 



(17a) 
(17b) 



18a) 



7?+(z) 

1 (E k -E ) 



where the self-energy S + (z) is give in Eq.flHJ). 

Substituting Eqs.(fT5|) into Eq.(fT7b). we then have 

1 



(18b) 



F(n) 



x 



7£ + T f 



i m T7 

AT ^ 

2 (Ek 



N ^ SI- E k + ie T]+(z)r]-(z) 



E^ 



0c 



0j i OT T 1 
+ ^J-dc-I-Oc- 



En 



(19) 



where rj~(z) is analytically continued from the lower half plane. As done in Eq.flBJ), trans- 
forming the summation over the wavenumber into the integral in the limit of N — > oo, the 
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explicit form of F(Q) is obtained as 

F(a) 



g 2 B 2 VB 2 - n 2 



(B 2 -Q 2 )(n-E ) 2 + g i B 4 



x 



(n-E 



Oc 



9 2 S 2 ■ ~^0c gB j • (20) 

Now we shall decompose the absorption spectrum F(Q) into the pole and the branch 
point contributions. For this purpose, we first rewrite F(Q) in terms of the contour integral 
by using the relation 



— / ,](,■■ ■■ 

2tt 



B 



dE k p(E k 



B 



1 



B 



— I dE k 

2m 



1 



B 



g 2 B 2 



[T 1 + {E k )-T 1 -{E k 



(21) 



where ■ ■ • represents a function of k. By taking the limit of N — > oo and applying these 
relations to Eq.([19jh F(Q) may be cast into the form of a contour integral: 



F(Q) = Im— / 
1 ' 2niJ T 



/ r il- E k + ier}(E k ) 

where the contour T is shown in Fig^a). As shown in FigJSJ the contour V can be deformed 
to that shown in Fig|3fb), where the cross denotes the pole location at E k = z d = E — ij. 
Along T, the solid and dashed lines are in the first and second Riemann sheets, respectively, 
so that r)(E k ) takes the corresponding value of r] + (E k ) and r]~(E k ) in Eq.()22|l. respectively 
In order to extract the pole contribution from Eq. (j22|) . we evaluate the residue around 
the pole, and obtain 



Fo(fi) 



Im — : / dE k 

2^ J pole 



1 



g 2 B 2 



X 



1 dc "r 1 Oc 



g 2 B 



n-E k + ier]+(E k ) 

2 , rp rp (0 - ^O) 

r ZirfcJoc" 



#5 



(23a) 



-Im 



iVrf 



, ^2 (^ _ -^o) 2 
J rfr "T J r 



Oc 



+2Trf c T ( 



Oc" 



where N d is given by 

n: 



d_ 

dz 
1 + 



(grf ~ E 
gB 

g 2 B 2 



En- 



z d {z d - Ep 
z 2 -B 2 



N ^ 

k 



E k 



-i z=z d 



(23b) 

(24a) 
(24b) 
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Subtracting the pole contribution from Eq. (j22j) . we can write the branch point contribution 

as 

FUQ) = Im— [ dE k - J— , (25) 

where the contour T' is depicted in FiglB^c). The total absorption spectrum is then decom- 
posed to F(p) = F (Q) + Fi(fi). 

Though the exponentially decaying unstable state corresponding to the pole in the second 
Riemann sheet cannot be identified in the Hilbert space, it is possible to identify it outside 
the Hilbert space. One of the authors (T.P.) et al. have shown that the Friedrichs solution 



can be decomposed into the unstable state 
(and (4>k\) in the generalized Hilbert space 



i) (and (4>d\) and the dressed field states \(j> k ) 



X>*><0*l = l&X&l+X>*>&l • (26) 

k k 

One can prove that this decomposition has one to one correspondence with the Fq and F\ 
spectral components. It should be noted that the factor Nd in Eqs. (l24|) is now recognized 
as a normalization constant of the unstable decaying state so that ((pdl^d) — 1, and then 
N d = (d\4> d )(<Pd\d). 



B. Fano Profile 



In order to clarify the enhancement of the branch point effect due to the singularity of 
the density of states in the absorption spectrum, we first show in Fig0]the calculated results 
for an artificial situation with T 0c = 0. We take B = 1.0 and T^ c = 1.0 in all calculations 
in the present work. In FigEJ the results for g = 0.2 are shown for (a) Eq = —0.1 and 
(b) E = —0.98. Solid lines are the total absorption F(Q) calculated by Eq.(j20J)? while the 
dashed and the chain lines are the pole contribution -Fo(^) i n Eq. (|23b ) and the branch point 
contribution Fi(Q) in Eq. (j23j) . respectively. In FigE^b), the spectra are magnified around 
Q ~ —B, while the overall spectrum is shown in the inset. 

When 1 — \E \/B ^> g (FigHfa)), the pole contribution F Q (Q) is dominant in F(Q); the 
Fi (Q) contribution is very small except for the tiny increase around the band edges. In this 
case, Nd — 1 and |i?o — -^ol gB in Eq. (j23b ). resulting in a sharp Lorentzian spectrum 
of F (Q) as shown in FigHfa). The fact that F (Q) is dominant in the entire energy region 
of —B < Q < B, suggests that the time evolution of the adatom localized state almost 
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completely obeys the exponential decay law, i.e. the Markovian process is predominant, as 
mentioned in the end of the previous section. 

As Eq gets close to the band edge (FigHJb)), the decay rate of the unstable state is non- 
analytically enhanced (as discussed in §2) and the energy shift of \E — E Q \ becomes large 
due to the singularity in the density of states of the miniband at the band edges, as shown 
in Fig. E| Therefore the spectral width of Fq(Q) gets wider and the shift of the peak position 
from Eq becomes visible as shown in FigEfb). Furthermore, it is found from Eq.([24b) that 
the divergence in the density of states enhances the second term in Eq. (|24b ). and thus Nd 
is largely reduced from 1 while Im[A^] remains negligibly small. Consequently, the relative 
ratio of the branch point contribution of Fi(Q) to the pole contribution F Q (Q) becomes 
larger than in FigEJ^a) . As discussed at the end of the previous section, the relative increase 
of the branch point contribution indicates that the time separation between the exponential 
and the power law decays becomes obscure in the time evolution of the surviving amplitude. 
Therefore, the non-Markovian process with memory effect becomes significant, which main- 
tains the quantum coherence in the decaying process. The reduction of Na (= (d\<f>d) (<f>d\d)) 
also indicates that the unstable state \<pd) (or (<pd\) contains a larger contribution from the 
electronic continuum components of \k). This suggests that the contribution of the contin- 
uum to the dressing cloud is more significant. As shown below, the persistence of quantum 
coherence in the decay process, or the large dressing effect, manifests the quantum interfer- 
ence between the different optical transition paths once the other absorption transition Tq c 
is introduced. 

Next we consider the actual case in which both optical transition paths are allowed: Tq c ^ 
and Tdc ^ 0. We show in Fig|3fa) the calculated results of F{VL) for the same parameters 
as in FigEfa) except that T$ c = 1.0 here. The pole contribution Fq(Q) and the branch point 
contribution F 1 (fl) are shown in Fig|5fb) and (c), respectively. All these spectra are depicted 
by the black solid lines. In each panel, the spectra are further decomposed into the spectral 
components due to the first, second, and the third terms in [• • •] of Eq. (f2T)j) or Eq.(J2Hb)- 
These terms are attributed to the d-d diagonal component, 0-0 diagonal component, and 
the interference term between these two in Eq. (|17b ). These are shown by the red, blue, and 
green lines, respectively. 

Introduction of T 0c changes the symmetric Lorentzian spectral shape of F(Q) shown in 
FigUJa) , and yields an asymmetric spectral shape around Q ~ Eq shown by the black line 
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in Fig|Sfa). As seen from Fig|3fa), this is caused by the interference term of Eq.()20|) shown 
by the green line in the figure. Around Q ~ ±B the F(Q) shows a sharp rise which reflects 
the divergence in the density of states, though F(Q) remains finite around Q ~ ±B and 
F(±B) = 0. It can be seen in FigJH^b) that F (Q) is the dominant contribution to F(Q) 
for Q ~ E and we see that the antisymmetric spectral shape of the interference component 
of Fq(Q) (green line) is the origin of the asymmetry in F(Q). In fact, neglecting Im[A^] 
in Eqs. fl23|) . which is appropriate except for \Eq\ > B, the interference term in Fq(Q) is 



approximately given by 



cyrp rp 

F inter {Q) = f^Oc^^ 

gB 

n — e e — e 



When the condition 1 — \E \/B 3> g is satisfied, it holds that Re[iVV] ~ 1 and E —E <^ 1. As 
a result, only the first term of Eq.()27p contributes to F ™ ter (FigJHfb), green line), leading to 
the antisymmetric spectral shape of the interference component of Fq(Q). The asymmetric 
spectral shape obtained when the branch point effect is neglected is represented by the 
well-known Beutler-Fano profile 2^|. 

However when Eq lies close to the band edges, \Eq — Eq\ becomes large, the second term 
of Eq. (|27b ) can no longer be neglected, which increases the interference component of F (Q) 
compared to the d-d diagonal component. We show in FigElthe calculated results for the 
same parameters as in FigEfb) where Eq = —0.98B, except that we take T<j c = 1.0. As seen 
in FigE^b) the interference component of F (Q) (green line) is enhanced compared with that 
in FigE^b). As mentioned above, in this case quantum coherence plays a key role in the 
decay process, which is clearly reflected through the enhancement of the interference effect 
in the absorption spectrum. 



IV. SUMMARY AND DISCUSSIONS 



In the present work, we have evaluated the charge transfer rate from an adatom impurity 
on a ID semiconductor superlattice. The decay rate is dramatically enhanced due to two 
square-root forms of singularities in the density of states. In the vicinity of the singularities 
at either edge of the band spectrum, the decay rate becomes a nonanalytic function of the 
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coupling constant g at g — as it can only be expanded in powers of g 4 / 3 . The time evolution 
of the localized adatom state is governed by the pole and branch point contributions to the 
decay rate, which account for the exponential and nonexponential decay, respectively. 

We have demonstrated that the absorption spectrum from the inner core level of the 
adatom to the conduction states is an appropriate probe to observe both contributions. 
While the branch point effect is usually small, it becomes important when the adatom lo- 
calized state is resonant with the miniband of the ID superlattice in the vicinity of the 
singularity in the density of states. In the case where the branch point effect is signifi- 
cant, the quantum coherence in the decay process becomes exaggerated. This explains the 
enhancement of the quantum interference in the absorption spectrum. 

The advantage of using the semiconductor superlattice is that it is easy to vary the pa- 
rameters in an electronic system. Using modern nanotechnology, we can vary the parameter 
values widely, in order to systematically investigate the effect on the decay process. Other 
spectroscopic techniques, such as the resonant optical light scattering spectrum, may also 
be used to investigate the time evolution of the unstable state in detail, which we are now 
studying. 
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FIG. 1: (a) An adatom attached to a ID quantum-dot array, and (b) level structures of the adatom 
localized state and the bound state in each quantum well. The width of each well is a few nm to 
100 nm. The adatom is located at the n = 0-th well. The two optical absorption transitions from 
the core level |c) with T^ c and Tq c are also shown by the arrows. 

FIG. 2: (a) Ei/B and Eq/B vs. Eq/B for g = 0.5, which are shown by the thick solid and dashed 
lines, respectively. The location of the critical values ±E^/B are indicated by the arrows. The 
thin line is y = Eq/B. (b) j/B vs. Eq/B for g = 0.1. The maximum value of ^y m ax/B occurs at 
Eq/B = ±1. 

FIG. 3: The contours of the integral for F(Q) Ea. (|22|) (a) and its deformation (b), and the integral 
for F (n) Eq.(|23i) ( c )- 

FIG. 4: The calculated F(Q) (solid line), F (Q) (dashed line), and Fi(Q) (chain line) for g = 0.2, 
and Tdc = 1-0 and Tq c = 0: (a) Eq = — 0.1-B and (b) Eq = —0.985. The thin vertical lines indicate 
the position of Eq. In (b), the horizontal axis is expanded around 17 ~ — B, while the overall 
spectrum is shown in the inset. 

FIG. 5: The calculated F(£l) (a), Fq(£1) (b), and iq(c) for the same parameters for FigOJa) except 
Tq c = 1.0: The spectra are decomposed into the d-d diagonal (red line), 0-0 diagonal (blue line), 
and the interference terms (green line). 

FIG. 6: The parameters here are the same as in FigOJb) except that Tq c = 1.0. 
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